function br = isParamsGood(pf, pStr2Double, Opt)
    global Opt; 
    
    t0 = 10000;
    ting = 365;                
    
    resParams0 = CalcResParams(pf, pStr2Double, 0);                   
    sol = NNsysExact_ing0(t0, resParams0);    
    v1 = deval(sol, t0);
    v2 = deval(Opt.optsol0, t0);
    d0 = abs(v1(1) - v2(1));
    dcomp = v2(1) * 0.3333;

    resParams15 = CalcResParams(pf, pStr2Double, 1.5e+6);        
    sol = NNsysExact_ingNonzero(ting, resParams15);    
    d15 = deviationlog10(sol, Opt.optsol15, ting);

    resParams75 = CalcResParams(pf, pStr2Double, 7.5e+6);
    sol = NNsysExact_ingNonzero(ting, resParams75);    
    d75 = deviationlog10(sol, Opt.optsol75, ting);

    br = [ (d0 < dcomp) (d15 < 0.27) (d75 < 0.27) ];  
end